General Localization Lengths for Two Interacting Particles in a Disordered Chain 
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The propagation of an interacting particle pair in a disordered chain is characterized by a set of 
localization lengths which we define. The localization lengths are computed by a new decimation 
algorithm and provide a more comprehensive picture of the two-particle propagation. We find that 
the interaction delocalizes predominantly the center-of-mass motion of the pair and use our approach 
to propose a consistent interpretation of the discrepancies between previous numerical results. 
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The problem of interacting electrons in a disordered 
potential is one of the important unsolved problems in 
condensed matter physics. This has been emphasized 
again by the recent observation Q] of a metal-insulator 
transition in two dimensional (2d) systems which was 
theoretically unanticipated. Some time ago, Shepelyan- 
sky Q proposed that it would be worthwhile to consider 
the simple case of two interacting particles in a random 
potential. He predicted that unexpectedly, such a par- 
ticle pair could propagate coherently over distances £2 
much larger than the single-particle localization length 
£1 as long as the two particles are within £1 from each 
other. 

Specifically, Shepelyansky obtained for the two- 
particle localization length 

6 ~ (u/w 2 ) 2 , (l) 

where U denotes the interaction strength and W the dis- 
order strength. Since £1 ~ 1/W 2 , Eq. (1) implies an 
enhancement of the localization length for weak disor- 
der. Shepelyansky 's original argument involved several 
uncontrolled assumptions for the single-particle eigen- 
states. This led to a number of (mostly numerical) at- 
tempts to study the problem of two interacting par- 
ticles more rigorously. Imry || rederived Shepelyansky's 
result, Eq. (1), by an extension of the Thouless block 
scaling picture. Frahm et al. j| computed £2 ~ W~ 3 3 
using the transfer matrix method (TMM). Von Oppen 
et al. § introduced a Green function approach, allowing 
one to project the problem on the subspace of doubly oc- 
cupied sites, and concluded £2 ~ U/W 4 . Subsequently, 
Song and Kim treated the idea of von Oppen et al. 
rigorously using the recursive Green function method and 
found £ 2 ~ W~ 2 9 . 

Recently, Romer and Schreiber [Q concluded from the 
TMM that the enhancement effect does not exist. In view 
of this claim and of the quantitatively different expres- 
sions for £2 quoted above, it appears that there are few 
secured results in this field. Our purpose in this paper 
is to present a more comprehensive picture of the two- 
particle propagation by defining and computing a set of 
localization lengths. We unambiguously show that the 
effect exists and propose a resolution of the controversy 
in the previous works (7ln(|. 



It is currently not clear whether these ideas have any 
relevance to the degenerate finite-density Fermi gas. It 
appears to be the most promising direction to consider 
the localization properties of quasiparticle pairs. There 
have been a number of studies ||]8j|j|] whether quasiparti- 
cle excitations delocalize relative to single-particle ones. 
While a numerical study for a one-dimensional (Id) sys- 
tem showed derealization only for unrcalistically high 
excitation energy of the pair (of order of the bandwidth) 
H , both arguments || and numerical studies || in higher 
dimensions suggest the possibility of a new pair mobility 
edge close to the ground state. 

The two-particle problem in one-dimension is de- 
scribed by the Hamiltonian 



H = H 1 ®l + l®H 1 + U^\m)\m){m\{T 



(2) 



where m labels the N sites of the Id lattice and H\ is 
the usual single-particle Anderson Hamiltonian 



[ e m \m)(m\ + t (\m)(m + 1| + \m+ l)(m|) 



(3) 



e m is a random site energy, drawn from a box distribution 
with — W/2 < e m < W/2, and U the on-site interaction. 
The hopping matrix element t is set to unity throughout 
this work. A convenient quantity to study the localiza- 
tion properties of the pair is the two-particle Green func- 
tion G = (E — 7Y) _1 . The two-particle localization length 
£2 on which previous studies have focused is defined in 
terms of G as || 

lim - — - — r < ln\(m,m\G\n,ri)\ >, (4) 
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where the double bracket denotes the disorder average. 

In this paper, we discuss general localization lengths 
which provide a much more comprehensive picture of the 
localization properties of the particle pair. First, we con- 
sider a general center-of-mass (CM) motion by defining 
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<C In I (m, 771 — a|G|77, n-a)| > . 

(5) 



1 



We find that, surprisingly, £2, a is essentially independent 
of the particle distance a, even if a exceeds the single- 
particle localization length £1. We also study the behav- 
ior of G for relative motion at fixed CM, as characterized 
by 

^jT 1 = — hm — <C In I (m + n, 7?i — n|G|m, to) I ^> . (6) 

n— >oo 71 

Finally, we consider the propagation of one of the parti- 
cles with the other held fixed, as described by 

£7* = — lim — <C In I (m, m + n\G\m, to) | 3> . (7) 
n — >oo 77 

As opposed to £2 and £2,0, we find that the latter two 
lengths are only very weakly affected by the interaction 
U. Nevertheless, it will turn out that these lengths are 
indispensable for obtaining a more comprehensive picture 
of the two-particle propagation and for understanding the 
discrepancies between previous numerical results. 

While £2 could be computed by projecting the problem 
on the subspace of doubly occupied sites, this is no longer 
possible for the generalized localization lengths defined 
above. For this reason, we introduce a new decimation 
algorithm, which allows us to compute these localization 
lengths efficiently. As opposed to the projection method 
for £2 used in ref. 0, this algorithm is numerically ex- 
act. We briefly describe the procedure for computing £2. 
Adaption to the other lengths defined above is straight- 
forward. Since the interaction acts only on symmetric 
states, we specify to (spinless) bosons. Using a sym- 
metrized basis 

. . ( \m)\m) if to = 77, . . 

\ mn) = \ (1/V2)(|to)|?7) + \n)\m)) xim^n. ^ > 

and interpreting (to, n) as sites of a 2d square lattice, 
the Hamiltonian of Eq. (0) can be interpreted as describ- 
ing a single particle on the 2d lattice shown by the thin 
solid lines in Fig. 1. The off-diagonal elements of 7i are 
nonzero only for nearest-neighbor bonds and equal to V2 
(1) if one (none) of the nearest-neighbor sites is a dou- 
bly occupied state. Our goal is to compute the Green 
function G(E) which is the inverse of a sparse matrix, 
D = E — H of linear size ~ TV 2 . Clearly, a direct manip- 
ulation of the whole matrix is inefficient both in terms of 
time and storage, and becomes forbidding for N > 100. 
To circumvent this problem, we recursively decimate the 
irrelevant matrix elements of the Green function. We 
start by decomposing Hilbcrt space into subspaces i, each 
of which is spanned by the states along one of the dashed 
lines in Fig. 1(a) and which are labeled by their dimen- 
sions 1 < i < N. We denote the projection of D onto 
these subspaces as Z>;. Clearly, D couples only neigh- 
boring subspaces (i) and (i + 1), and we call the cor- 
responding (7 x (i + l) dimensional) coupling block in 

(n) 

the Hamiltonian V%. Finally, we define vectors x^ with 
elements 



(*r)i={N-i+j,j\G{E)\nn), (9) 

given by matrix elements of the Green function G be- 
tween a doubly occupied site \nn) and the states in sub- 
space i. Since only neighboring subspaces are coupled, 
one readily derives from DG — 1 the set of coupled lin- 
ear equations 

D,4 l) + = 0, 

Ci x l-i + Ax,'" 1 + =0, 2 < i < N - 1, 

V%-iy$-i + D n xl%> = en, (10) 

where e„ is the N dimensional unit vector with (e„) m = 
&n,m- Solving these equations, we obtain 

XiV = SNGn, (11) 

where the Qi can be computed recursively from 

Qi = (A ~ V^Gi-iVi-x)- 1 with Qi = Di 1 . (12) 

Finally noting that (x^ 1 ); = (ll\G(E)\nn), we can now 
compute the localization length £2 from Eq. (Q). This 
reduces the calculation to manipulations of matrices of 
sizes from 1 x 1 to N x N. It is worthwhile to point out 
that at the final stage of the iteration, the calculation is 
formally reduced to an effective Id model for a single par- 
ticle. It is straightforward to generalize the algorithm to 
compute the other localization lengths defined in this pa- 
per. E.g., £ r is calculated by decomposing Hilbert space 
according to the dotted lines of Fig. 1(a). 

For £2) we set n = 1 in the above algorithm and obtain 
= <C In |(Zi|G|ll)| > with 1 < I < N for each param- 
eter set (W,U). We find that t;,i depends linearly on I, 
implying an exponential decay of the Green function. To 
eliminate finite size effects near I = 1 and I — N, we fit 
t/,i in the range N/5 < I < 4iV/5 to 

ti,i = —r+c (13) 

with c a constant. We find that for chains N > 200, 
our results for £2 are essentially independent of system 
size N, suggesting that finite-size effects on £2 are rather 
weak. Similar procedures are performed for the other 
localization lengths. 

Our main results are presented in Fig. 2. All data have 
been obtained for system size iV > 200 and for the cen- 
ter of the band E = 0. In view of the special nature of 
the doubly-occupied sites due to the on-site interaction, 
it is natural to ask whether the definition for the two- 
particle localization length £2 correctly captures the CM 
motion. To answer this question, we plot £2,0 for inter- 
action strength U = 1.0 as function of a in Fig. 2(a). We 
find that £2,0 remains unchanged up to rather large a, 
implying that £ 2 is indeed a good description of the CM 
motion. In fact, £2,0 remains independent of a even for 
a > £1. Hence, as opposed to the previous beliefs [|]||] 
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the interaction affects the two-particle motion even if the 
particle distance exceeds £1. This can be understood in 
terms of single-particle propagation in the 2d lattice of 
Fig. 1. £2, a is associated with the transition probability 
along the dashed line a distance ~ a from the diagonal. 
We recall that when (to, to — a\G\n,n — a) is expanded in 
powers of the hopping matrix element t, it is given by a 
sum over all possible paths from \n, n — a) to \m, m — a). 
If the distance between these two sites, ~ |m — n\, is 
much smaller than a, the effect of the interaction would 
be negligible. However, £ 2l a is defined by the limiting 
behavior of |m — n\ — > 00 with a finite, cf., Eq. (5). In 
this case, the contributions of paths which are sensitive 
to the interaction U are no longer negligible and £2,0 re- 
mains influenced by the interaction even though a > £1. 

In Fig. 2(b), we show our results for £ r (symbols) which 
describes the decay of the Green function with relative 
distance. For comparison, we also plot £2 (lines). At 
[7 = 0, the two lengths are equal within the numerical 
accuracy, i.e., £ 2 — £r- As U increases, £ r remains nearly 
constant while £2 shows a pronounced enhancement in 
qualitative agreement with Shepelyansky's prediction Q . 

In Fig. 2(c) we plot £/ (symbols) and £ 2 (lines), where 
the former describes the range over which one particle 
moves with the other one fixed. At U = 0, we find that 
£f approximately equals to 2£ 2 - As already seen in Fig. 
2(b), £ 2 shows a strong increase with U. By comparison, 
£/ shows a much weaker increase. Hence, there exists a 
U C (W) beyond which £2 exceeds £/. 

At U = 0, the two particles move independently and 
the propagation of a given particle is not affected by 
whether the other is moving in the reverse (£ r ) or in the 
same direction (£ 2 ), implying £ r — £ 2 . Moreover, since £2 
measures two-particle propagation, while the single- 
particle motion, one expects that the transition proba- 
bility for £2 is given by the square of that for so that 
£2 = We note in passing that &(U — 0) 7^ £i/2, as 

was pointed out in ref. ||. Since both £ r and £/ are de- 
termined by the limiting behavior for diverging distance 
between the particles, one does not expect them to be 
strongly influenced by the interaction. This explains the 
rather flat dependences of both lengths on U. 

With the additional information from £ r and £/, we 
can now construct a wavefunction picture in the 2d lat- 
tice representation of the problem (Fig. 1). At U = 0, 
this implies that the wavefunction profile is described by 
a square as shown by a thick solid line in Fig. 1(a). As 
U increases, the length of the edge associated with £2 in- 
creases while that associated with £ r remains essentially 
constant. For U > U c , the wavefunction profile becomes 
highly anisotropic and we find that it can be well de- 
scribed by an ellipse as shown by the thick solid line in 
Fig. 1(b). The elliptical shape predicts the relation 

= 26k 

We have checked that our data are in good agreement 



with this expression for U > U c . This clearly shows that 
the enhancement effect is associated predominantly with 
the direction of £2, i-e., the CM motion of the two parti- 
cles. 

These results allow us to resolve some of the above 
mentioned discrepancies between previous numerical 
studies. We start by noting that the TMM measures 
the largest length scale from the iV 2 Green-function en- 
tries (ln\G\Nm) with 1 < to, n < N fTc| ], According 
to our results, there are two competing lengths £2 and 
£/. For U < U c , we find that the largest length is 
while for U > U c , it is £ 2 - Therefore, the TMM actu- 
ally measures £/ for U < U c and £2 only for U > U c . 
We first compare our results to those of Frahm et al. [Q . 
We find that their results for the localization length are 
two to three times larger than our result for max{£/, £2} 
at given values of W and U . We attribute this to large 
finite-size effects in the TMM, as suggested in ref. 0. On 
the other hand, in ref. 0, any enhancement effect was 
attributed to finite-size effects and it was suggested that 
the TMM produces the single-particle localization length 
for a sufficiently large system size. This is clearly incon- 
sistent with the results of the present paper. We suggest 
the following explanation for the numerical results of ref. 
@ . The argument of ref. @ is based on TMM data for 
(W,U) = (3.0,1.0). For these parameters, we find that 
the largest length is £f = 13.2 ± 0.3, which is close to 
£1 ~ 11.7. Hence, we expect that the data in ref. 
in fact extrapolate to £/ which is indistinguishable from 
£1 within the numerical accuracy of ref. W. Therefore, 
we contend that the principal argument of ref. is a 
misinterpretation of data for a special parameter set and 
expect that the TMM exhibits the enhancement effect 
clearly once £2 3> £/■ Finally, we find that the ^2's in 
ref. [j5| are somewhat larger than those in this paper, 
which we attribute to the approximate treatment of the 
Green function in ref. Iq]. 

We thank H.A. Weidenmiiller for useful discussions. 
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(b) U > U c 

FIG. 1. Sketch of the two dimensional lattice (thin solid 
lines) and the wavefunction profile (thick solid lines). The 
dashed lines (dotted lines) represent the index scheme for the 
calculation of £2, a (£r). Lengths are measured in terms of the 
lattice constant d. The factor \/2 arises because £ r and £2 are 
defined in units of the diagonal length of the smallest square 
of the lattice (V2d) while is defined in units of its edge 
length (d). 
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FIG. 2. (a) £2,0 as function of a. (b) £2 (lines) and £ r 
(symbols) as function of U. (c) £2 (lines) and £/ (symbols) 
as function of U. All data have been obtained for system size 
200 < N < 300 and for E = 0. 
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